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' Abstract. We consider a stochastic functional differential equation with an arbitrary Lipschitz 

diffusion coefficient depending on the past. The drift part contains a term with superlinear 
$_j ' growth and satisfying a dissipativity condition. We prove tightness and Feller property of the 

segment process to show existence of an invariant measure. 
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1. Introduction and preliminaries 



There have been quite some investigations on stationary solutions of stochastic functional dif- 
ferential equations with nonlinear diffusion coefficients, see for instance [1, 3, 9] and references 
therein. One approach is to rewrite the functional differential equation as a semilinear infinite 
dimensional equation and use results on invariant measures of such equations (see [5]). The 
operator induced by the linear part of a functional differential equation is often not dissipative. 
For results on invariant measures for non-dissipative systems, see [2, 12]. These results require 
that the linear part is exponentially stable and that the Lipschitz constant of the diffusion is 
small with respect to the decay of the linear part. By means of a finite dimensional analysis it 
, has been shown that the Lipschitz constant of the diffusion coefficient may be arbitrary large, 

provided the diffusion coefficient is uniformly bounded (see [8]). 
<T> 1 In this paper we prove existence of an invariant measure for stochastic functional differential 

equations with no boundedness conditions on the diffusion coefficient nor conditions on the size 
of its Lipschitz constant. Instead, we consider a stabilizing feedback term in the drift with 
superlinear growth. Let r > and denote by C([— r, 0],R ) the space of R d valued continuous 
functions on [-r,0] and let g: C([-r,0],R d ) -» R d and h: C([-r, 0], R d ) R dxm be Lipschitz 
functions with respect to the maximum norm. Let (B(t))t>o denote a standard R m -valued 
Brownian motion defined on a filtered probability space (f2, (^)t,P). We will show existence 
of an invariant measure for the functional differential equation 



dx(t) = [ -x(t) ■ \x(t)\ s + g(x t ))dt + h(x t )dB(t), t>0, (1.1) 

where s > and xt denotes the segment of x given by 

x t (9) = x(t + 9), 0e[-r,O]. 

In order to show existence of an invariant measure, we consider the segments of a solution. In 
contrast to the scalar solution process, the process of segments is a Markov process. We show 
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that the process of segments is also Feller and that there exists a solution of which the segments 
are tight. Then we apply the Krylov-Bogoliubov method. 

Since the segment process has values in the infinite dimensional space C([-r, 0},R d ), bound- 
edness in probability does not automatically imply tightness. For solution processes of infinite 
dimensional equations, one often uses compactness of the orbits of the underlying determin- 
istic equation to obtain tightness. For an infinite dimensional formulation of the functional 
differential equation, however, such a compactness property does not hold. 
Our proof of tightness involves a Lyapunov function technique to obtain boundedness in proba- 
bility for the segment process (xt)oo- Further we use the assumption on the coefficients for the 
deterministic part, and Kolmogorov's criterion for the noise part. By using a monotonicity argu- 
ment we prove the Feller property for (xt)t>o which implies the existence of an invariant measure 
by the Krylov-Bogoliubov Theorem. Our analysis holds true for the more general equation 

f dx(t) = (f(x(t)) + g(xtj)dt + h(x t )dB(t), for t > 0, 
[ x(s) = tp(s) for s E [— r, 0], 

where we assume the following hypotheses: 
(H ) / : R d -> R d is continuous and 

hm — j — — = — oo. 

H-H-oo \v\ z 



(Hi) g : C([-r,Q],R d ) -» R d , h : C([-r,0],R d ) -> R dxm are continuous and bounded on 

bounded subsets of C([-r, 0], R d ). 
(H2) There exists a positive constant L such that for all x, y G C([— r, 0],] 



(2(/(x(0)) - /(y(0)), x(0) - y(0))+ + 2(g(x) - g(y), x(0) - y(0)) + 

+ \lh(x)-h(y)\f)<L\\x-y\\ 2 , 

where |M| := (Tr(MM*)) x / 2 denotes the trace norm of the matrix M. 

The initial process ip has almost surely continuous paths and is independent of (B(t)) t >o with 
E\\^(-,uj)\\p < 00 for aUp> 2. 

Note that under hypotheses (Ho), (Hi) and (H2) and thanks to [10, Theorem 2.3], equation 
(1.2) has a unique global solution given by 

x(t) = x(0)+ [ f(x(s))ds+ [ g(x s )ds+ [ h(x s ) dB(s) for any t > 0. 
Jo Jo Jo 

We will prove existence of an invariant measure fi for the segment process (xt)t>o associated 

to the solution x(t)t>o- Of course our hypotheses (Hi) and (H2) allow the coefficient h to be 

degenerate which can not guarantee uniqueness of \i. For recent results on the uniqueness of 

invariant measures for stochastic functional differential equations, see [6]. 

We end this introduction by the following elementary remark which is useful for our arguments 
in the sequel of this paper. 

Remark 1.1. Let T > 0. Consider a stochastic process x(t), —r<t<T with continuous paths 
and let xt, t > be its associated segment process on [— r, 0]. If xq = ip and p>l, then 

E sup \\x t \\ p < E\\(p\\ p + E sup \x(t)\ p 

0<t<T 0<t<T 
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Proof. We have 

E sup ||xt|| p = E sup sup \x(t + s)\ p 

0<t<T 0<t<T -r<s<0 

= E sup sup |x(s) | P 

0<t<T t-r<s<t 

= E sup \x(s)\ p <E\\(p\\ p + E sup \x(s)\ p . 

-r<s<T 0<s<T 



2. Tightness of the segment process (x t ) t > 

In this section we will prove tightness of the family {xt : t > 0}. To this end we shall prove first 
boundedness in probability. 

We fix the initial process ip and consider the solution of (1.2). 

Proposition 2.1. Under hypotheses (Ho), (Hi) and (H2) the process (xt)t>o is bounded in 
probability. 

For the proof of the proposition we need some preparation. Let n: [0,oo) x Q — > R be a 
progressively measurable process with locally square integrable sample paths. Consider a one- 
dimensional Brownian motion (/3(i))t>o and for fj, > let us introduce the following equation 

J dv(t) = -nv{t)dt + rj(t, u)d/3(t), t > 
I v(0) = 0. 

If we denote by (f^(-)) its solution we have 

v,{t)= f e-^- s \{s,u)d(3{s). 
Jo 

The following lemma gives an estimate for the process w M (-)- 

Lemma 2.2. For 2 < p < +00 and fi > 0, there exists a positive constant a p ^ such that 

lim a p . = 

and 

E sup \vn(t)\ p < a, p ^ -E f \t](s,uj)\ p ds, for every T > 0. (2.1) 

0<t<T Jo 

Proof. Fix 2 < p < 00, T > and assume that E J Q T \t](s,uj)\ p ds < 00. Let | < a < \ and 
define 

y(t):= [ (t-s)' a e'^ t ~^7 1 (s,Lo)d(3(s), t > 0. 
Jo 

Using the factorization formula (see [4, Sect. 7.1]) 

f e-rt-»)r,(s,u) dp(s) = S -^^R a y(t) 
Jo n 

where 

R a f(t)= [\t-sr- 1 e-^f(s)ds 
Jo 
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defines a bounded linear operator from L p ([0,T],R) into C([0, T], R). Indeed, take a function / 
in D>([0,T],R), then we have 

\R a f(t)\< /V-sr-v^i/wid* 

JO 

n-i 

ft 



< 



< 



Therefore 



E sup 

\0<t<T 



f 

Jo 



3 -/i(t-s) 



Lf([0,T],]E 
Lf([0,T],I 
LP([0,T],f 

77(5,0;) 



/ (*-*) 

JO 

f + OO 

Jo 



(a-l)p/(p-l) e -w(t-s)/(p-l) ds 



p — 1 



l)p/(p-l) e -//ps/(p-l) ^ s 



ap — 1 



p-i 



p — 1 



E sup 

\0<t<T 



sin 7ra 



7T 



< \\R a \\ (E||y^"" 



LP([0,T], 



< 



p — 1 
VP 



r 



ap — 1 
p — 1 



E||y(-)lli P([0 ,T], ffi ))"- 



Using Burkholder-Davis-Gundy's inequality we obtain 



Elll/ll? 



LP([0,T],R) 



;/ |y(t)|^i 
Jo 

= / E| f\t-s)- a e-^ t -^r t (s,uj)dp(s)\ p dt 
Jo Jo 

< c P E^ T QT(t - s)- 2a |e-^- s )r ? (s,w)| 2 ds^) ' dt 

rp E rp 

(Young's inequality) < c p Qf s - 2a e _2/M ds^* • E J \v(s,uj)\ p ds 



Hence we have 



milt 



LP([0,T],R) ^ C P ■ ( v (2jU) l- 2a r ( 1 " 2 «)J E _/ W a > w )l P ds 



c p,m • E / 1^0, ^)I P ds, 



where c PiM := c p ( ( 2/J )i-^ T(l — 2a)^ 2 . Therefore we deduce 



E 



sup / 

\0<t<T JO 



3 -/x(t-s) 



77(5,0;) d/3(s) 



where 



P-l 



pa— 1 / 1 \ p— 1 

7 — 1 \ 



p-l 
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We are now in the position to complete the proof of the proposition. 
Proof. Let A > 1. For x € W 1 we define 

R x (x) :=2(f(x),x)+X\x\ 2 . 
By hypothesis (H ) there exists A\ > such that 

(f(x),x) 
\x\ z 

Again by (Ho) we can find B > independent of A such that 

R\(x)<B + \A 2 x forallx€M d . (2.2) 

We now consider the solution x(-) of equation (1.2) and set z(t) := |x(i)| 2 , t > 0. Then Ito's 
formula implies that for fixed t > we have 

dz(t) = 2(f(x(t)),x(t))dt + 2{g(x t ),x(t))dt + \lh(x t )fdt + 2(x(t), h(x t )dB(t)) 

= ( - Az(t) + R x (x(t)) + 2< 5 (x t ), x(i))cft + \\h(x t )fyt + 2(x(i), h(x t )dB(t)) 

< ( - A*(t) + fl A (x(t)) + 2(g(x t ) - 5 (0), s(t)) + 2(g(0),x(t)} 

+ 2|/»(x t ) - h(0)\f + 2\\h(0) I 2 ) dt + 2(x(t), h(x t )dB(t)) 

< ( - Xz(t) + R x (x(t)) + 2L||x t || 2 + 2(g(0), x(t)) + 2||/i(0)|| 2 )di + 2{x(t), h{x t )dB{t)) 

< ( - Xz(t) + fl A (x(t)) + 3L||x t || 2 + i| 5 (0)| 2 + 2||/i(0)|| 2 )di + 2(x(i), h(x t )dB(t)), 

(2.3) 

where we used the estimate 

(g(0),x(t)) < ^\x(t)\ 2 + ^| 5 (0)| 2 < ^||x t || 2 + ^b(o)| 2 . 

Set D := ^^(O)] 2 + 2||/i(0)|| 2 , so the variation of constants formula yields 
z(t) < z{0)e~ xt + e- x(t -^ (r x {x(s)) + 3L||x s || 2 + D) ds + 2 J* e ~ x(t ~ s ^ (x(s), h(x s )dB{s)) 
< z ( ) e ^t + Al + ^±R + EL sup \ x (s)\ 2 + 2 f e~ x(t ' s \x(s),h(x s )dB(s)). 

A A - r <s<t JO 

There exists a one-dimensional Brownian motion (5 with respect to the same filtration such that 

(x(s),h(x s )dB(s)) = rj(s,uj) d(3(s), 

where 

j=l i=l 

By (H 2 ), we get 

|r,(^)| 3 < \x(s)\ 3 \jh(x s )\f < 4 |x( S )| 3 (L 3 / 2 ||x s || 3 + D 3 / 2 ). (2.4) 
Hence for < t < r we obtain 

e t z(t)<z{0)+e r (Al + ^-^-)+^e r sup |x(s)| 2 + 2e r sup f e^^^s) d(3(s 

^ A / A - r <s<t 0<t<r JO 
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Now using Lemma 2.2 and (2.4) we get 



E sup 

0<t<r 



e- x ^ V (s)dp(s) <a 3 , x rE\\ Vr 

o 



3 



< 



|3 



4ayr (D^EII^f + L^^EHxrf+Edl^fllxrll 3 ))) 
< 2a 3iA r(Z) 3/2 (E||x r || 6 + l) + L 3/2 (3E||x r f +E\\(p\ 



If we choose k € (1, e 3r ) and 7 > 1 such that (a + b + c + d) 3 < kcl 3 + 7(6 3 + c 3 + d 3 ) for all 
a, b, c, d > we have 



E sup \e l z{t)\ 3 < kE\z(0)\ 3 + -/e 3r (A 2 x + —— ) + 7— ^— e 3r (E||(/9|| 6 + E sup \e s z(s)\ 3 ) 

0<t<r ^ AS 0<s<r 

+ 7l6a 3iA re 3r ( J D 3/2 (E sup |e'z(t)| 3 + l) + L 3/2 (3E sup |e'z(i)| 3 + E||y?f) Y 

^ 0<t<r 0<t<r ' 

Let ip{s) := \(p(s)\ 2 , s G [-r,0]. We define the function F : C([-r,0],R) -» M+ by 



F(C):= sup (e 3s \((s)\ 3 ). 

-r<s<0 

We deduce from the above calculation that 

B + D\3 27L 3 



(2.5) 



EV(z r ) < K e- 3r EV(^j) + 7 (a\ + — ±— ) + 7^e 3r (EF (Y>) + EV(z r )) 

+ 16 7 a 3j A r (EV(z r )e 3r {D 3 ' 2 + 3 L 3 / 2 ) + e 3r L 3 l 2 EV(^) + D 3 ' 2 ) . 
Hence, for A* sufficiently large, we get 

EV{z r ) < SEVOj) + p, (2.6) 



where 

K e~ 3r + 7^#e 3r + I67 a 3 \ re 3r L 3 ' 2 
1 - 7e 3r (^f + 16 a 3 ,A. r (D 3 / 2 + 3L^)) ' 

9 '~ 1 - 7 e 3 '- + 16 o 3 ,a. r (DV 2 + 3L 3 / 2 )) ' 

provided that EU(z r ) < 00 (EV(ip) is finite by assumption). To see that this property holds, 
apply the previous calculation to the process \x(i)\ stopped as soon as it reaches level N and 
then let N — > 00. Iterating (2.6) we get 

EV{z kr ) < 5 k EV(ip) + — —c < EV(tp) + — —c , for all fc € N. (2.7) 

l—o l—o 

Let i > 0. Then there exists k € No such that kr < t < (k + l)r and we have 

E|M| 3 <E||z fcr || 3 + E||z (fc+1)r || 3 . (2.8) 

Using (2.7) we obtain 

E||z fcr || 3 = E sup \z kr {s)\ 3 < e 3r EV{z kr ) < e 3r (EV(ip) + 9 

Combining this with (2.8) yields 



1-5 



supE||x t || 6 < +00. (2.9) 

t>o 
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This implies the boundedness in probability of the segment process (xt)t>o and the proposition 
is proved. ■ 

The following theorem is our main result in this section. 

Theorem 2.3. Under hypotheses (Ho), (Hi), (H2) the family {C(xt), t > 0} is tight. 

Proof. From (2.9) we have in particular the boundedness in probability of the finite dimensional 
process (x(t))t>o and hence the family {C(x(t)), t > 0} is tight. To prove the theorem, it is 
sufficient to show that 

/ \ 

= for any 7 > 0. 



lim sup P 



sup \x(v) — x(u)\ > 7 

t<u<v<t+r 



(2.10) 



To shorten notation let 



Thus we can write 



9(V) ■= 9(V) + f(rj(0)), r/€C([-r,0],] 



x(t)=x(0)+ [ g(x s )ds+ f h(x s )dB(s) 
Jo Jo 



and we have 



sup \x(v) — x(u)\ > 7 < P sup / \g( 

>,<v<t + r j \ t<u<v<t + r J u 

v~u<S 



t<u<v<t+ 

y v— u<8 



r )\ds>^ 



I 



+ 

= M t + N t . 



sup 

t<u<v<t+r 



rv 

/ h(x s ) dB(t 

J u 



> 



1 



(2.11) 



Let e, R > 0. For the term M t we have 

/ 



M t < 



sup / \g(x s )\ds>- 

<u<v<t + r J u ^ 



\ 



t<u<v<t + r J u 

y v— u<s 



\xt\\ < R, ||^t+r-|| < R 



J 



+ F({\\x t \\> R}) +F({\\x t+r \\> R}). 

Since the process (xt)t>o is bounded in probability we can choose R so large such that 

P({|H| > i?}) + P({||x t+r || > R}) < e - for all t > 0. 

By (Hi), g(x s ), s G [t — r, t + r] is bounded on the set {||xt|| < R} fl {||xt +r || < R}, so it follows 
that there exists Sq > such that 

( r - 

P sup / |5(x s )| ds > — \\xt\\ < R, \\xt+ r \\ < R I = for any 5 < Sq. 

t<u<v<t + r J u ^ 

y v— u<5 



Therefore we get 



lim sup M t = 0. 

S^o t >o 
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For the term N t we define 

J(t) := [ h(x s )dB(s). 
Jo 

Using Burkholder's inequality and (H2), we get 

E\J(t) - J(s)f =E| / h{x u )dB{u 



< cE 



\\h{x u )fdu 



(Jensen's inequality) < c\t — s| 3 (supE||x u || 6 + l), 

where c depends on L and D. Using (2.9) and Kolmogorov's tightness criterion (see [7, 2.4.11] 
or [11]) we infer that 



lim sup Nt = lim sup 1 

<5->o *>o 5_>0 t>o 



sup 

t<u<v<t-{-r 

y v— u<s 



h(x s ) dB(s] 



> 



7 



= 0. 



This establishes (2.10) and the proof is complete. 

3. INVARIANT MEASURES 



In this section we discuss the existence of an invariant measure fi for the segment process (xt)t>o- 
Since in the last section we proved tightness of this process, in order to apply Krylov-Bogoliubov's 
theorem we need to prove the Feller property of (xt)t>o- 

Proposition 3.1. Assume hypotheses (H ), (Hi) and (H 2 ). Let (i^ m ) me N be a sequence in 



C([— r, 0],R ) such that ip m ► ip. Let x m (resp. x) be the solutions to (1.2) with initial 

condition ip m (resp. if). Then for any t > 0, 



', sup \x m (s) — x(s)\ ^0 as m — > +00. 

t-r<s<t 



(3.1) 



(3.2) 



In particular, (xt)t>o is a Feller process. 

Proof. Using Ito's formula we can write 

d\x m (t) - x(t)\ 2 = 2(f(x m (t)) - f(x(t)) + g(x?) - g(x t ),x m (t) - x(t))dt 
+ \jh(x?)-h(x t )\fdt + dM(t), 
where 

M{t) := [ 2(x m (s)-x(s),(h(xT)-h{x s ))dB(s)} 
Jo 

is a martingale with quadratic variation process bounded by 4L J Q * \\x™ — x s \\ A ds. Thus if we 
define M*(t) = supM(s) we obtain 



s<t 



This implies 



IK m - x t \\ 2 < \\<p m - ip\\ 2 + L [ \\xf - x s \\ 2 ds + M*(t). 

Jo 

E||xr - ^|| 4 < 3E (\\<p m - iff + L 2 (£ \\ x ? ~ Xsf + (M*(t))^ 

<?>(\W m -¥\\ A + L 2 t E||x™ -x s \\ A ds + AL E||x™ - x t \\ A ds 



(3.3) 
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Hence, by Gronwall's inequality, we obtain 

Hx? - x t t < 3\\<p m - ^|| 4 e 12L '+ 3L2 ' 2 . 
This implies in particular that for tp : C([— r, 0], M. d ) — > M. d bounded and continuous we have 

lim Etp((x m )t) = E,ip(x t ) for any t > 0, 

which yields the Feller property. ■ 

Now, by the Krylov-Bogoliubov Theorem (see Sect. 3.1 in [5]) we have the following result. 

Theorem 3.2. Under hypotheses (Ho), (Hi) and (H2) the segment process (xt)t>o correspond- 
ing to (1.2) has an invariant measure. 

Remark 3.3. Our proofs show that hypothesis (Ho) can be weakened by requiring that 

lim sup ^ppf^ < —A, for A sufficiently large positive constant (which depends on L, r and h(0) ). 
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